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ABSTRACT 

We analyse the radial and tangential velocity fields in the Galaxy as seen from the Sun by using as a first approximation a simple 
axisymmetric model, which we then compare with the corresponding fields in a barred Af-body model of the Milky Way. This provides 
a global description of these quantities missing in the literature, showing where they take large values susceptible to be used in future 
observations even for sources well beyond the Galactic center. Absolute largest proper motions occur at a distance slightly behind the 
Galactic Center, which are there 1.5 times larger than the highest local proper motions due to the Galactic differential rotation. Large 
proper motions well beyond the Galactic center are well within the current astrometric accuracy. 
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1. Introduction 

In view of the ongoing and forthcoming advances of astrom- 
etry applied to the Milky Way well beyond the Solar neigh- 
borhood, such as the measurements of the trigonometric par- 
allaxes and proper motions of distant stellar masers (e.g., Reid 
et al. 12009b , or the upcoming GAIA survey, we develop here 
a simple first-order model of a galactic disk making explicit 
where large proper motions and radial velocities due to differen- 
tial Galactic rotation can be expected in the Galaxy. It is indeed 
important to develop quantitative insight on these observables 
and separate the effects due to differential rotation from those 
due to, for example, a bar. 

The level of analysis is elementary and could figure in as- 
tronomy textbooks. Despite literature search and discussions 
with astrometrists the points raised here do not seem to be 
[ commonly known. The probable reason for the lack of a pre- 
vious similar discussion is that observing proper motions has 
been considered for a long time as difficult and limited to a 
few kpc. Further, dust extinction has also limited observations 
in the Galaxy plane in the optical. Therefore, a local analysis 
has been considered as sufficient (e.g. Mihalas [19681 Mihalas & 
BinneyHSSH Binney & Merrmeld[I998]). 

In radio astronomy, however, the need to develop a global 
vision of the Galaxy has been reached much sooner. An analysis 
close to the present one concerning the Galaxy radial velocity 
field has been already done (see e.g. Burton [T9741) , since extinc- 
tion in neutral hydrogen is weak and integrated radial velocity 
field can be directly measured over most of the Galaxy disk. 

The most interesting aspect of the present analysis concerns 
the proper motion field at distances well beyond the Solar neigh- 
borhood, which shows a discontinuous behavior just beyond the 
Galactic center allowing to measure easily stellar proper mo- 
tions there and well beyond. The reason is easy to understand 
if we consider that the disk beyond the Galactic center is mov- 
ing in first approximation with a velocity opposite to the Sun 
velocity. The linear velocity difference of about 2 x 220 km s" 1 
largely exceeds the inversely to distance decreasing proper mo- 
tions, leading to proper motions up to about 440 km s _1 /9 kpc « 
10 mas yr" 1 , a well measurable quantity with current techniques. 



Of course the circular velocity is not equal to the average 
azimuthal velocity, the velocity that can be observed by aver- 
aging proper motions, due to the non-negligible contribution of 
the velocity dispersion to radial support. Especially close to the 
Galactic center the velocity dispersion increases to values com- 
parable to 220/ V3kms _1 « 130 km s" 1 , therefore one should 
see, as regions beyond the Galactic bulge are measured at in- 
creasing distances, a progressive change from individual random 
motions at about 130 km s _1 /8 kpc « 3 mas yr -1 superposed to a 
systematic proper motion with respect to an inertial frame due 
to the Sun motion of about 220 km s _1 /8 kpc « 5 mas yr -1 , pro- 
gressively increasing to about 400 km s -1 / 10 kpc « 8masyr _1 
as we explore regions 2 kpc beyond the Galactic center. This in- 
sight was already exposed by Binney (1995) in a little read paper 
concerned primarily with the non-axisymmetric motion effects 
induced by a bar. He already showed that large proper motions 
can be expected near the Galactic center even in the absence of 
a bar. In any case, the systematics introduced by differential mo- 
tion in the bar/bulge region must not be omitted, as pointed out 
by Zhao et al. (119951 . 

We also check that a fully consistent model of the Milky Way 
including the effects of a bar and velocity dispersions does not 
change the main conclusions reached with the simple model. For 
this, we use one of the Milky Way Af-body models of Fux (11997b 
including a bar and spiral arms to quantify the differences exist- 
ing between a fully consistent, self-gravitating, stable and barred 
model of the Milky Way, and the simple axisymmetric model. 

The plan of the paper is as follows. In Section 2 we give some 
definitions. In Section 3 we develop the first order axisymmetric 
model with analytic tools and a more elaborated rotation curve 
model. In Section 4 we analyze Fux's A/-body model and com- 
pare it with the axisymmetric ones and with some observational 
results. In Section 5 we present our conclusions. 



2. Decomposition of a differential velocity field 

Contrary to the tradition, in the following we express formulas 
as much as possible in Cartesian coordinates in order to avoid 
trigonometric expressions masking simple geometrical interpre- 
tation. Bold symbols denote vectors. 
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Fig. 1. The differential velocity field for h = 1, p = -0.05, 
components. 



*o - 8, yo = (left) and its decomposition in radial (middle) and tangential (right) 



For any velocity field V(x) e R 3 , x e R 3 one can define the 
differential velocity field relative to the velocity of a particular 
observer located at xq, moving at the velocity of the field V(xq), 



v(x) = V(x)-V(x ). 



(1) 



The scalar radial velocity field v r (x) with respect to the observer, 
the usually observable quantity derived from the Doppler shift, 
is the projection of v(x) along the line of sight unit vector 



n r (x) = (x - xo)/d, 



(2) 



where d = |x - xol is the distance to the observer location. The 
vector velocity field v r (x) is proportional to v r and directed in 
the radial direction n r . 



v r (x) = v(x) • n r , 



v r (x) = v r (x) n r . 



(3) 



The tangential velocity field is the difference between the 
velocity field and the radial velocity field, 



V/(x) = v(x) - v r (x). 



(4) 



The observable proper motion field ji is the tangential veloc- 
ity field V/(x) divided by the distance d to the observer 



Vr(x) 



(5) 



This is a 3 -dimensional vector field orthogonal to the line-of- 
sight, tangential to the celestial sphere. It can be converted to two 
angular velocity components for the desired spherical coordinate 
system using standard spherical trigonometry formulas. Here we 
will only use the magnitude of yu. 



3. Axisymmetric models 

The purpose of an axisymmetric model is to provide some in- 
sight into rotation curves resembling the Galaxy one. Therefore 
we choose first a particular representative but simple model of 
the Milky Way rotation curve that can be handled analytically. 
In a second step we use a more realistic but still axisymmetric 
model of the Milky Way (Kalberla[2003]). 



3.1. Simple model 

We consider a simple extension of the popular logarithmic ax- 
isymmetric galactic potential. We use cylindrical coordinates 
(R,z) (R 2 = x 2 + y 2 , with the origin at the Galactic center 
x = y = z = 0). The potential has the form: 



!(>+fH) 

|ln(l + S + |) 



for p*0, 
for /; = 0, 



(6) 



where p is a parameter, p e [-1,2], allowing to study different 
rotation curves going from an asymptotic Keplerian model (p = 
-1), to a constant rotation curve (p — > 0) (the usual logarithmic 
potential), to an harmonic potential (p = +2). The constants h 
and zo denote the respective horizontal and vertical scale lengths. 
The circular velocity reads 



WW 



1/2 



R 



1 Rl 
1 + F 



p-2 



(7) 



When p - the rotation curve v c (R) is asymptotically flat 
(for R » h), while for negative/positive values of p it de- 
creases/increases. For p = 2, it reduces to a straight line, a 
solid rotation corresponding to a constant density mass model. 
For p = - 1 it decreases asymptotically like a Keplerian model 
where the mass is concentrated at the center. When, additionally 
to p = 0, we send h —> 0, the rotation curves becomes a con- 
stant v c = vo for all R > 0, which corresponds to Mestel's disk 
rotation curve. 

We could have chosen to start directly with a rotation curve 
expression which represents the effective average azimuthal ro- 
tation field. However not any rotation field is consistent with a 
solution of Jeans' equations and a positive mass distribution. The 
advantage of the circular rotation curve derived from a potential 
is that it gives a rotation field consistent with a specific potential- 
density pair, while its drawback is that the circular rotation speed 
overestimates the effective mean rotation speed. 

The circular velocity vector field for a clockwise circular ro- 
tation like in the Milky Way reads, 



V(x) = -v c n z An R = {y, -x, 0} , 

K 



(8) 



where n# = {x,y, 0}/R, and n z = {0, 0, 1} are unit vectors in the 
galactic plane and perpendicular to it, respectively. The differ- 
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Fig. 3. On the left, contour map of the magnitude of the proper motion \i as observed from the observer's position at xo = 8, yo = as in Fig.Q] 
The dark/light grays correspond to low/large values. On the right, the side view of the same contour map shows the relative importance of proper 
motion, with the most negative, largest absolute value slightly beyond the Galactic center as seen from the observer's position. 



ential vector fields v(x), v r (x), and \ t (x), as well as //(x), can be 
derived in analytic form using the definitions given above. 
The radial vector field reads 



v r (x) xy - yx 



vo 



hd 



1 R1 

1 + ¥ 



h 2 



(9) 



where d is the distance to the Sun. This expression vanishes 
along the line xyo = yxo, and along the circle R = Ro. 
The tangential vector field reads 



v,(x) xxo+yyo-R 2 



vo 



hd 



, R 
1+ ¥ 



- + ~^ 1+ | M10) 



Here the location of vanishing transverse velocity is not simple. 
Vf(x) is symmetric with respect to the line xyo = yxo. At the 



limit d — > the field is proportional to d, that is, continuous but 
not difFerentiable. This is a well know result from the classical 
textbook analysis of the local motions around the Sun, where 
in term of Oort's constants A = \{v c /R - dv c /dR)\R and B = 



-j(v c /R + dv c /dR)\ Ro , we have 
v t = d(Acos(2/) + £). 



(11) 



This quantity converges to at d — > 0, but with a slope depend- 
ing on /. 

For the sake of illustration we show some figures for relevant 
cases. We set h = 1 and p = -0.05 to study conditions similar 
to those occurring in our Galaxy for a population of objects with 
a positive velocity dispersion, which sees its effective rotation 
decreased by the asymmetric drift. We put the observer at the 
position xo = {8, 0, 0}, which, in kpc unit, corresponds approxi- 
mately to the Sun position with respect to the Galactic center. 
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Fig. 4. On the left, rotation curve for the analytic model (dashed red line, same parameters as in Fig. [TJ and the Kalberla model (solid blue line). 
On the right, proper motion fi along the line passing through the Galactic center as observed from the observer's position at R = 8, in mas yr -1 
for v = 220 km s _1 and kpc length unit, for the analytic model (dashed red line, same parameters as in Fig.Q]) and the Kalberla model (solid blue 
line). 
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The differential velocity field is shown in the left panel of 
Fig. [T] and the radial component is shown in the middle panel 
of Fig. [T]and in the left panel of [2] Along the x-axis, where the 
observer is placed, the radial velocity vanishes, as well as along 
the circular orbit of the observer at radius R 8, as noticed 
above. The maximum radial velocities within the observer orbit 
are located on both sides of the Galactic center direction, and at 
distances slightly shorter than the Galaxy center. This field has 
been used extensively in radio astronomy where the neutral hy- 
drogen radial velocity can be measured over most of the Galaxy. 

The tangential velocity field is shown in the right panels of 
Figs.[T]and[2] As can be expected the tangential velocity is small 
in the direction of the Galaxy center on the same side as the ob- 
server, because the matter turns in the same direction at a similar 
speed there, while the side of the Galaxy beyond the Galactic 
center sees large tangential velocities because the matter turns 
in opposite directions. But since only the proper motions are ob- 
servable, we must consider the proper motion field, shown in 
Figs. El 

There is a strongly negative proper motion at the point 
(x = -1.395, y = 0) in a region behind the Galactic center, a 
region characterized by large amplitude absolute proper motion, 
fi ^ -8.4 mas yr" 1 for vo = 220 kms" 1 and kpc length unit, 1.5 
times larger than the highest local proper motion contribution 
due to galactic rotation (ji -5.6 mas yr -1 ). For positive values 
of the x coordinate, the proper motion ji is nearly zero, except 
at the observer position where a sharp discontinuity occurs. This 
follows from the above remark about the local tangential differ- 
ential velocity close to the Sun (Eq. (11)). Since [i = v t /d, then 
lim^o J~i = A cos(2/) + B, i.e., is not single-valued. 

The horizontal proper motion as a function of distance along 
the line Sun-Galactic center takes the following form, 



yU (x) = -■ 



vo 



1 + 



h 2 



p_2 



Xo 

~h 



v»2 \ 4 



1 + 



h 2 



(12) 



The zeros and extrema of no(x) do not have simple analytic ex- 
pressions. This curve is shown in Fig. [4] (right panel, dashed red 
line) for a representative choice of the free parameters. Clearly 
the proper motion maximum amplitude occurs for objects lo- 
cated slightly beyond the Galactic center. 



3.2. Kalberla's model 

Here we check with a more elaborated rotation curve model that 
our basic result holds, i.e. that the proper motions slightly be- 
yond the Galactic center have largest amplitudes. 



We take the model of Kalberla (12003b and extract 40 data 
points from the rotation curve figure. These data are then repre- 
sented as a cubic B- spline curve in the symbolic software pack- 
age Maple 13. This B-spline object can be considered as a usual 
function on which standard operations, including differentiation, 
can be applied. The analysis can be repeated exactly as for the 
previous model. 

Since the results turn out to be qualitatively the same, we 
only show the initial rotation curve and the final proper motion 
curve along the direction of the Galactic center in Figs. 0] (solid 
blue lines). The main difference occurs near the center inside 
the inner 3 kpc where the rotation curve presents a bump and 
a steeper raise. This feature may be due to the streaming mo- 
tions in the bar and, for applying it to observed stars, it should 
be corrected from the large velocity dispersion there by decreas- 
ing the average circular motion. Thus the higher proper motion 
amplitude (up to 12 mas yr -1 ) obtained here beyond the Galactic 
center is probably an upper bound of what could be observed. 

The differences between the models are largest precisely 
in the bulge-bar region where the absolute proper motions are 
largest, i.e., easier to obtain. Proper motions therefore represent 
a powerful observable for discriminating models of the bulge- 
bar kinematics. 
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Fig. 5. Density distribution for the barred Af-body model m08t3200. 
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Fig. 6. Rotation curve for the barred Af-body model m08t3200 (solid 
blue line) and the analytic model (dashed red line, same parameters as 
in Fig.Q}. 



4. Barred TV-body model 

The previous models based on rotation curves do not take into 
account the decrease of average rotation due to asymmetric drift, 
nor the effects of a bar. Both effects are largest in the bulge-bar 
regions, therefore it is useful to compare our previous models to 
a fully self-consistent model of the Milky Way including a bar. 

Here we analyse one of the barred models of the Milky Way 
obtained by Fux (11997b and we compare it with the axisymmet- 
ric models discussed in the previous section and with observa- 
tional results. 

Fux's barred models result from the self-consistent Af-body 
evolution of a bar-unstable axisymmetric model, with 10 5 parti- 
cles in the nucleus- spheroid component (which is well suited to 
represent the nuclear bulge and the stellar halo) and in the disk 
component, and 10 5 particles in the dark halo component (which 
has been added to ensure a flat rotation curve at large radii). 
The spatial location of the observer has been constrained by 
the K-bznd map, obtained by the Diffuse Infrared Background 
Experiment (DIRBE) on board the COBE satellite and corrected 
for extinction by dust, while masses have been scaled accord- 
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Fig. 7. Contour map of the magnitude of the differential velocity field 
for the barred Af-body model m08t3200 decomposed in its radial (top) 
and tangential (bottom) components. 



ing to the observed radial velocity dispersion of M giants in 
Baade's window. Among the models obtained by imposing such 
constraints, the model named m®8t320® is recommended in 
Fux (1997) since it reproduces satisfactorily the kinematics of 
disk and halo stars in the Solar neighborhood. 

The density distribution of m®8t320® is shown in Fig. [51 By 
a rotation of the particles the observer has been moved to the 
position xo = {8, 0, 0} in kpc unit. The angle between the line 
joining the observer to the Galactic center and the major axis 
of the bar equals to = 25° and the corotation radius is Rl = 
4.8 kpc (Fux [T997l 

The rotation curve is shown in Fig. [6l solid blue line. For 
comparison, we also show the rotation curve for the analytic 
model discussed in the previous section, for p = -0.05, h = 1. 
As in Kalberla's model, the main difference occurs near the cen- 
tral 3 kpc, with the absence of a bump. 

The differential velocity field v defined in Eq. (Q]) is now ob- 
tained by averaging the particle velocities V* located at different 
z positions to smooth out irregularities due to granularity, and by 
subtracting the averaged velocity vector V(xo) of the particles 
in a volume of 1 kpc 3 centered at the position of the observer 
xo. The radial velocity field and the tangential velocity field are 
shown in Fig.[7](top and bottom panels, respectively) and should 
be compared with the fields shown in Fig. [2] for the analytic ax- 
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Fig. 8. Contour map of the magnitude of the proper motion \i as observed from the observer's position at x = 8, y = (left) and proper motion 
jio along the line passing through the Galactic center (right) for the barred Af-body model m08t3200. 
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Fig. 9. Distributions of radial velocities (left) and proper motions (middle) obtained by selecting stars with < / < 20° and positive/negative x in 
the barred Af-body model m08t3200. Right panel: proper-motion distributions obtained by selecting stars with < / < 20° and different distance 
d from the observer. 



isymmetric model. From these panels it can be seen that, apart 
from noise in the boundary cells due to the finite number of par- 
ticles, the overall structure of the contour maps is in agreement 
with those corresponding to the analytic model. The presence 
of the bar does not introduce asymmetries changing the overall 
results. In order to quantify the effect of the bar, we have con- 
structed a symmetric system by randomising the particle polar 
angles. Thus, we could evaluate that the relative error in both ra- 
dial and tangential velocity fields due to the presence of the bar 
is less than 15%. 

The proper motion field (averaged on z) and the proper mo- 
tion yuo (averaged on z and \y\ < 1 kpc) are shown in Fig. [8] 
(left and right panels, respectively). Also in this case the over- 
all structure of the contour map agrees with the analytic model 
map (compare with Figs. [3]) and the behavior of jiq is analog 
to that shown in Fig. [4| (right panel) apart from numerical os- 
cillations due to the presence of the singular point at the ob- 
server's position. The relative error between the barred model 
and the one obtained by randomising angles is 6%. As in the ax- 
isymmetric models, in Fux's barred model a region of particles 
with large-amplitude proper motion exists, located behind the 
Galactic center. In this region, the proper motion is of the order 

of// 7 masyr" 1 for particles at distances to the observer of 

10 kpc. The same order of magnitude holds in regions slightly 
above and below the Galactic plane, which can be reached by di- 
rect measurements. This order of magnitude for the proper mo- 



tion is well within the astrometric accuracy of the GAIA survey, 
which will be 12-25 //as at the 15 th magnitude and 100 - 300 //as 
at the 20 th magnitude (see e.g. Brown 20Q7])- 

We compare now the predictions of this A/-body model with 
some observational results in order to understand how the pa- 
rameters of the Galactic bar can be constrained by measurements 
of radial velocities and proper motions. Mao & Paczynski (2002 ) 
used two simple analytic models to show that samples of stars in 
the far and near sides of the bar display differences in both radial 
and tangential streaming motions. In that paper, it was suggested 
to use red clump giants to constrain the kinematics of stars in 
the Galactic bar. Red clump giants are bulge stars which occupy 
a distinct region in the color-magnitude diagram. They have a 
well-defined peak in their observed luminosity function, which 
can be used to select stars at the bright slope of the peak, and thus 
closer to us, from those at the faint slope which are more dis- 
tant. Following this suggestion, data of the Optical Gravitational 
Lensing Experiment II (OGLE-II) in the Baade's window were 
used to derive the difference in the average proper motion be- 
tween the bright and faint samples of red clump giants (Sumi et 
al. 12003b and the result A(//> -1.5 masyr -1 is consistent with 
the value estimated in Mao & Paczynski (12002b who assumed 
a streaming motion of the bar of v ~ 100 kms" 1 in the same 
direction as the solar rotation. 

More recently, different methods were used in order to se- 
lect bulge samples. Vieira et al. (2007) measured proper mo- 
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tions inside Plaut's window with an accuracy of 1 masyr -1 us- 
ing ground-based data. They selected a bulge sample at the 
mean distance of 6.37 kpc in the red giant branch by cross- 
referencing with the Two Micron All Sky Survey catalog. 
Clarkson et al. (12008b reported on the use of the Advanced 
Camera for Surveys WFC on the Hubble Space Telescope to 
extract proper motions in the Sagittarius window with an accu- 
racy of 0.3 mas yr" 1 . They selected kinematically a bulge sample 
by introducing a cutoff on longitudinal proper motions and on 
proper-motion measurement errors, as suggested by Kuijken & 
Rich ( 2002). An interesting result of such studies is that longi- 
tudinal proper motions can clearly separate near and far samples 
of stars. 

We use Fux's accurate model in order to evaluate the effect 
of distance on radial- velocity and proper-motion measurements. 
We consider almost 24 '000 particles with Galactic longitude 
< / < 20° and we compare their distribution in the near region 
x > with that in the far region x < 0, averaged on \z\ < 0.5 kpc. 
The distributions of radial velocities v r and proper motions [i 
are shown in the left and middle panels of Fig. [9l respectively. 
From the left panel, it can be seen that the far sample in x < is 
slightly shifted toward more negative radial velocities, in qual- 
itatively agreement with Mao & Paczynski (2002). The shift 
between the two samples is A(v r ) ~ 5.4 kms -1 . Longitudinal 
proper motions show much clearer association with distance 
than radial velocities. The middle panel in Fig. [9] shows that the 
distribution of far- sample proper motions is shifted toward more 
negative values with respect to that of the near sample, with a 
shift of A(/i) -2.8 mas yr -1 . Moreover, the dispersion of the far 
sample is 60% smaller than that of the near sample. 

In order to discuss the effects of contamination of bulge sam- 
ples by disk stars we select stars in Fux's model with different 
mean distances from the observer. Since the corotation radius is 
Rl = 4.8 kpc and axis ratio is b/a = 0.5 in Fux's model, sam- 
ples in < / < 20° and 6 < d < 10 (in kpc length unit) consist 
mainly of bulge stars, while those in d < 6 or d > 10 have 
large contamination by disk stars. In the right panel of Fig. [9] we 
show the proper-motion distributions obtained by selecting stars 
within these regions. By comparing the right and middle panels 
of Fig. [9l we see that the effect of the contamination from dis- 
tant ( d > 10) and nearby (d < 6) disk stars is to move the mean 
value of the bulge samples toward more positive values. Distant 
disk stars also have the effect of reducing the dispersion of bulge 
proper-motion distribution, as discussed in Vieira et al. (120071) . 
Knowledge of the global distribution of proper motions as shown 
in Figs. [3] and [8] help us to understand and to evaluate these ef- 
fects of contamination. 

5. Conclusions 

Starting from the analysis of simple axisymmetric models of the 
Milky Way, we have studied the differential velocity field, and 
in particular the proper motion field, providing a global descrip- 
tion of observables. Interestingly, there is a region behind the 
Galactic center extending to several kpc further out where the 
proper motions induced by differential rotation have the largest 
amplitudes in all the Milky Way, an apparently little known basic 
feature (with the exception of Binney 1995) with direct observa- 
tional consequences. Here we clarify on a global scale the effects 
related to first order to the average circular rotation and to second 
order to a bar. 

We have checked that these high proper motions are also 
found in self-consistent barred configurations, such as the Milky 
Way Af-body model of Fux (119971) . The proper motions at dis- 



tances of 10 kpc to the observer are of the order of 7 masyr -1 
and thus are well within the accuracy of the forthcoming surveys 
such as GAIA. It is crucial to be able to separate the sources ac- 
cording to distance, at least to separate the sources closer to or 
beyond the Galactic center. 

Provided that the position of the sources, such as stars or 
OH masers, are intense enough to be measured with respect to 
an absolute frame of reference, as well as their radial velocities, 
the global velocity field of the Milky Way including non-circular 
motions should become accessible. This will be a first order in- 
formation for constraining the strength of the bar and the mass 
distribution in the bulge, the keystone of the whole Galaxy. 
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